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FIRST  LINE  OF  TEXT 


THERMOACOUSTIC AL  INSTABILITY  OF  A  HETEROGENEOUS  GAS  FLOW 

FIRST  Llt.E  OF  TITLE  | 

K.  I.  Artamenov  and^ I.  G.  Krutikova 


The  Investigation  of  the  stability  of  a  one-dimensional  single- 

i  '  * 

velocity  flow  with  distributed  internal  heat  sources  reduces  to  finding 

the  characteristic  values  of  the  boundary- layer  problem  for  a  system  J 

of  three  ordinary,  differential  equations.  j 

i  i 

The  characteristic  values  have  been  numerically  determined  on 

j 

digital  electronic  computers  for  two  laws  of  heat  conduction:  l)  the 

i  ■  I 

heat  supply  rate  is  proportional  to  the  density  of  the  gas,  and  2)*  the; 

i  1 

heat  supply  rate  is  proportional  to  the  mass  of  the  gas  in  a  given  j 

cross  section  of  the  flow.  J  J 

For  a  flow  with  a  temperature  gradient  varying  along  its  length  | 
(first  case)  the  interaction  of  the  pressure  and  entropy  waves  leads  J 

to  the  excitation  of  inherent  acoustical  oscillations  in  the  flow;  the 

!  ! 

fundamental  note  is  the  first  to  become  unstable.  The  stability  ! 

i  I 

boundary  is  determined  by  the  values  of  three  dimensionless  criteria:  | 
the-ratio  of  the  specific  heats  7,  the  Mach  number  at  the  inlet  Mer;-antl 

the  criterion  p,  characterizing  the '  relationship  between  the  heal  ' 

• - 2  '  !  2 - ! 

supplied  to  the  gas  and  the  initialjheat  content  of  the  gas.  Th<f  lastj 

driter!on-deterjttipe£.£both--the  -  steady-state-  temper  ^distribution-  —  -I 
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along  the  length  and  the  relationship  between  the  heat  supply  rate 

r  L!Ni  CF  7c  A  7  ; 

I  t 

and  the  density  of  the  gas  In  the  presence  of  oscillations. 

A  flow  with  constant  temperature  gradient  along  its  length  (second 

1  .  i  ,  : 

case)  proves  stable  in  the  investigated  intervals  of  t>\e  criteria.  i 

I  ! 

Autoexcitation  of  oscillation  in  gas  flows  is  a  fairly  frequent 

i  !  i 

occurrence,  an  example  of  which  could  be  the  acoustical  instability  j 

! 

of  the  operational  process  in  the  combustion  chambers  of  Jet  engines,  j 
The  mechanism  exciting  the  oscillations  consists  in  the  interdependence 


c  :*T  ■  •  jP  C  f  T  Tl  r 


between  the  rate  at  which  mass  or  energy  is  supplied  to  the  flow  and  j 

the  velocity  and  thermodynamic  parameters  of  the  flow.  This  inter-  j 

l 

dependence  is  usually  realized  in  certain  zones,  the  dimensions  of  1 

i 

\ 

which  are  small  in  comparison  to  the  dimensions  of  the  flow  as  a  whole; 

i  1 

and,  consequently,  to  the  wave  length  of  the  oscillations. 

t  | 

In  the  schematization  of  the  problem  it  is  possible  to  introduce  j 
certain  surfaces  with  well  defined  characteristics  which  separate  the  ^ 

!  I 

flow  into  regions  for  each  of  which  the  equations  of  homogeneous  flow 

i  i 

are  applicable.  When  a  one-dimensional  flow  is  examined  in  this  way,  i 

the  pressure  and  entropy  disturbances  generated  by  the  dividing  surfaces 
are  propagated  Independently  of  one  another  and  may  interact  only  at  J 
the  boundaries  of  the  homogeneous  regions.  | 

Under  certain  conditions  these 1  oscillations  may  increase  through  I 

:  !  i 

time.  The  problems  of  flow  stability  for  this  schematization  have  been 

I  -  -  : 

examined  at  length  in  Rauschenbach1 s  book  [1],  • 

■  i  ! 

In  a  number  of  cases  this  schematization  is  unsatisfactory  and  t 

;  ;  i 

it- is- necessary  to  take  into  account  the  continuous  variation  of5the — ! 

parameters  of  the  main  flow.  An  example  might  be  tfie  ^ 

I - 3  !  3 - j 

problem  of  the  stability  of  a  one- dimensional  flow  with  the  heat2 _ j 

source  i  distributed  along  its  entire  j  length.  For  a  description  of  the  -j 
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oscillations 'in  such  a  flow,  it  is  necessary  to  make  use  of  the  com-  ; 

f-T.'.r  Cr 

plete  energy  equation,  which  expresses  the  continuous  interdependence! 
between  the  pressure  waves  propagating  with  the  velocity  of  sound  ! 

■  .  i 

(acoustical  waves)  and  entropy  waves  propagating  at  a  velocity  equal  i 
to  the  flow  velocity.  We  note  that  such  a  flow  may  be  considered  j 

concentrated  when  col/u  «  1,  but  not  when  co i/c  «  1  as  in  the  i&oen- 
tropic  case  (u  and  £  are  the  velocity  of  the  flow  and  the  velocity  of 

I 

sound  respectively;  <d  is  the  oscillation  frequency;  and  L  is  the  ; 

characteristic  length  of r the  f  low)  ^  j 

I  I 

1.  Equations  of  Heterogeneous  Plow.  Consider  a  one-dimensional, ‘ 

flow  of  an  ideal  gas  in  a  channel  with  varying  cross-sectional  area,  i 
Internal  heat  sources  are  present  in  the  gas;  the  heat  supply  rate  ' 

i 

per  unit  length  is  denoted  by  Q.  -I 

!  i 

When  customary  notations  are  used  the  flow  equations  have  the 

following  form:  ' 

•  '  ! 

_  -  ■  •  A 

dm 


1  W  ,  9Fm  - 
7-5T  +  -3T  =  °. 


_ dz  t-i n 

**  9*  *  dt  Tpdt  ~F~  W 


(1.1) 


The  solution  to  the  equations  will  be  sought  in  the  form 

!  I 

*  =  />#(l+»l).  P  =  P°{l+«),  u  =  b°(1+v),  ^  =  (1.2) 

' 

t 

The  quantities  with  the  superscript  *  characterizing  the  funda- 

! 

mental  steady-state  flow  regime  are  found  from  Eq.  (l.l),  if  we  dis- 

i  t 

card  the  first  term  on  the  right  side  in  the  operator 

'  i 

|  1 

d  9  ,  9 

!  dt  “  dt  ‘  *  a. 


dt 


-Hie  dimensionless  small  disturbances  r\,6,v,  and  £  satisfy  aA Bystem 


! 


of  three  linear  partial  differential  equations: 

*+■•£+£)-<>.  £+--S+£2+£(*+*-«-* 


-l 

JL_ 


2-T%  +  *-&-rS)-a<9-r,-r)  .(B-j 
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Here  ft  is  the  characteristic  frequency  of  the  heat  supply  process. 

~f  :$7 \  K.S  C? 

When  ft  «  0  the  entropy  equatioh  is  solved  separately  from  the  other 

I 

two;  if  SI  ji  0,  then  this  quantity  determines  the  degree  to  which  the 

I 

entropy  oscillations  are  interdependent  upon  the  other  flow  parameters. 

I  i 

We  note  that  for  isoentropic  distributed  mass  and  momentum  trans¬ 
fer  the  flow  would  be  described  by  a  system  of  second  order  differen- 

| 

tial  equations,  instead  of  third  order  as  in  our  case. 

We  shall  limit  our  consideration  to  a  flow  with  small  acceleration; 
i.e.,  we  shall  s6t 'duV^z'  —  'O.  In  this  case  the  area  of  the  channel  ! 

i  ! 

cross  section  is  determined  by  the  heat  supply  law.  Let  us  consider  j 

the  two  heat  transfer  laws:  j 

1  i 

the  heat  supply  rate  is  proportional  to  the  density  of  the  gas  J 

!  -  ! 

C  variant  l)  (1.4)  j 

i 

I  ! 

the  heat  supply  rate  is  proportional  to  the  mass  of  gas  in  a 

I  ; 

given  cross  section  i  1 

I  ' 

Q  =  ktpF  ( variant  2)  (1.5) 

1  I 

I 

We  shall  assume  that  the  static  heat  supply  law  is  preserved  in  I 

the  presence  of  oscillations;  then^  we  will  obtain  for  both  cases:  j 

i  I 


(l.6> 


The  coefficients  in  the  equations  are  not  explicit  functions  of 
time,  and  therefore  the  solution  may  be  sought  in  the  form 


6  =  6  (*)  r*  5 

_Let  us  introduce  the  dimensionless  coordinate  £  *  z/L  and  t^e. 


dimensionless,  time  t  =  tuo/L;  then  Eqs.  (1.3)  may  be  reduced  to2 


r 


>rd 


canonical  form 


1- 

0. 
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Here 


(T  _  j »/«)".  = « (_  r  _  7x)  +  v  (—  a/v  +  rM  +  n  (rA#*Xr  +  r*) 
(T  -  M')  -3|-  -  W  +  V  (A/V  -  T*)  +  1  (—  r  —  T*) 

(T  —  A/*)  =  —  *rA/*X  +  v  (—  rA/V  +  T-V*X)  +  n  (A/V  +  r^A/*) 


(1.7) 


A/*  =  — ?»— 
P*/P*  * 


X  =  Q  -  =  -  ^  , 

s  p 


The  quantity  M  in  our  case  is  the  isothemic  Mach  number. 

The  coefficients  M*  and  X  are  determined  by  the  heat  supply  law. 
Making  use  of  the  equation  of  state  of  an  ideal  gas  we  obtain  from 
the  steady- state  equations  j 


Vi  +  vl 


_ 

2(i+-pl) 


2  QJ> 

r  -  p.u.f.ct. 


(1.8) 


for  variant  1,  and 


P^«Cp  T, 


(1.9) 


for  variant  2. 

!  ! 

We  use  the  subscript  0  to  denote  the  values  at  the  inlet  (when 

i 

£  m  o) .  For  the  boundary  conditions  we  take  the  following: 


•V  +  6  ='  0,  tj  —  76  =  0  whmn  g  =  0 

v  +  Vi*  -  V*  *1  =  0  ,*«»  £  =  i 


(1.10) 


The  first  and  second  conditions  express  the  constancy  of  flow  rat^ 
and  entropy  respectively  at  the  iniet  to  the  channel.  The  third  con- i 

I  i 

dition  expresses  the  constancy  of  the  Mach  number  at  the  outlet  and 

I  1 

corresponds  to  quasi3teady  discharge  through  a  critical  nozzle.  j 

The  boundary  condition  at  the  inlet  into  a  supersonic  nozzl^  was~^ 

| - A  I  4  I 

given .more  precisely  by  Crocco  [2]^  we  shall  confine  ourselves  tp  the | 
simple  relationships  listed  above,  which  are  valid  when  the  length - | 

of  the  convergent  section  of  the  nozzle  is  small  in  comparison  to  the  ^ 

l _ a _ i _ o _ j 
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over-all  length  of  the  channel. 

t  f ; «  .  ~-r  Ct  1 :  A. . 

2.  Method  of  Solution. 


'  T ' 


I  .  ! 

The  problem  of  Investigating  the  flow  i 


stability  reduces  Itself  to  a  determination  of  the  condition8  for 

i  ! 

which  the  real  part  of  the  characteristic  values  of  the  boundary- layer 
problem  (l.7)-(l.l0)  pats  through  zero.  When  Re  r  >  0  the  flow  is  j 
unstable^  when  Re  r  <  0  it  is  stable  .  ‘  j 

i  ; 

The  characteristic  values  of  the  problem  in  seme  domain  D  of  the  1 

I  ■ 

complex  variable  r  -  x  +  iy  aref oundin_  the_  following  way.  The  | 

f  I  *^r  re  Tl 'r‘  r  l  I 

domain  D  is  covered  with  a  network  of  values  of  x  and  y.  For  each  of1 

I  ! 

these  values  Cauchy’s  problem  is  solved  for  syBtem  (1.7)  with  the  j 

I  1 

initial  conditions  for  €  -  0.  i  ! 


v  6  =  0,  tj  —  76  =  0,  6  =  1 


(2.1) 


Naturally,  with  these  initial  conditions  the  third  boundary 

I 

condition  in  ( 1 . 10)  was  not  fulfilled;  when  £  =  1  we  obtain: 


v  -r  V*  ■ 6  —  7*  1  =  ti  (x,  y)  -f  Ci  (x.  y) 


(2.2) 


The  characteristic  numbers  of  the  boundary-value  problem  under  | 
consideration  are  those  values  of  r  *  x  +  iy  for  which  simultaneously^ 


£1  (x,  y)  =  0,  c*  (x,  y)  =  0 


(2.5) 


The  calculation  of  Cauchy's  problem  was  programmed  on  a  Strela  1 

,  1  1 

electronic  computer  (subsequently  other  machines  were  used  also).  j 

I  I 

The  characteristic  values  were  j  determined  graphically  for  chosen  j 

values  of  7,  Mo,  and  M-  using  the  values  of  C*(x,y)  and  Cz(x,y)  Obtained 
on  the  machine.  The  accuracy  of  the  characteristic  values  was  4e^  ^ 

‘ - 3  I  3 - 1 

^termined  by  the  accuracy  of  the  solution  to  the  Cauchy  problem.  2  The_j 

solution  to  system  (1.7)  Is  fairly ;  complicated  in  form.  The  furictloni 

L _ A _ 1 - 0 _ j 
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6(4)  is  a  rapidly  oscillating  function  since  a  term  containing  e“  ^ 

fif;:  r:  cf  ”>•  ,  j 

is  included  in  its  solution  (the  period  of  the  oscillations  is  2ir/y)., 

The  function  v(4)  also  varies  rapidly,  but  its  shape  is  smoother  ! 

'  ! 

than  that  of  6(  4)  since  the  coefficient  of  v(  4)  in  the  second  equation 

2  '  { 

of  (1.7)  is  Mo  times  less  than  that  of  6(4)  in  the  first  equation. 

The  function  q( 4)  varies  little  by  comparison  to  6( 4)  and  v( 4) »  since 

St 

the  coefficients  in  the  equation  for  tj(4)  are  approximately  Mo  times: 

1  I 

less  than  the  coefficients  in  the  equations  for  6( 4)  and  v( 4)  .  In  j 

addition  to  these,'  ’a  numberrof  the  coefficients  in  the  equations  are  I 

1  | 

far  from  smooth  (the  sharp  maximum  when  4*0  for  X  when  p.  is  large,  j 

..  | 

etc!) .  1  1 

The  equations  were  integrated  over  each  division  interval  of  the  J 

1  >l 

4  axis  in  the  following  way.  1 

The  function  6(4)  was  determined  by  exact  integration  of  the 
first  equation  in  (2.7)  with  q( 4)  replaced  by  a  constant  and  v( 4)  I 

j 

replaced  by  a  linear  function.  The  function  v(  4)  was  determined  from 


the  second  equation  by  exact  integration  after  substitution  of  6( 4)  J 

found  previously.  The  equation  for  tj(  4)  was  computed  using  '  ; 

!  ! 

i  4 

Fuler's  method  on  an  average  layers  1 


The  method  taken  for  the  calculation  possesses  slow  convergence  , 
for  large  and  therefore  for  acceptable  accuracy  the  number  of  | 

i  | 

division  intervals  N  for  the  majority  of  variants  was  chosen  in 


-  2 .»/„  ,nL(*  + 


y  =  An  - 


1  —  u  * 


M 


T  —  1  \ 

(3.1) 

"  2T  ) 

0.  1. . .  .) 

(3.2) 

When  M*  «  1  the  frequency  of  the  oscillations  corresponds  tp  the 
characteristic  acoustical  frequencies  of  the  tube: 

a  =  M_v  =  =  fa, 


Krc. 


(3.3) 


/ 
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Fig.  1. 


For  p.  =  0  but  when  the  second  condition  is  changed  to  tj  —  7,6'=  0, 
7i  ^  7  for  £  =  0,  there  occur  in  addition  to  the  ” fundamental"  values 
of  (3.2),  a  large  number  of  characteristic  numbers  for  which 
y«!2mr(n  =1,2...). 

There  is  also  an  analogous  exact  solution  for  small  p.  when  the 
coefficients  of  Eqs.  (2.7)  are  replaced  by  certain  average  values. 
Obviously,  the  new  characteristic  values  having  periods  with  respect 
to  £  equal  to  2nrr  correspond  to  the  entropy  waves.  As  evaluations 
have  shown,  the  damping  coefficients  for  them  for  the  values  of  p.  that 
have  been  studied  lie  well  into  the  negative  region  (x^  —  3),  there¬ 
fore  these  characteristic  values  were  not  taken  into  account  when 
finding  the  stability  regions. 
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It 'follows  from  Eqs.  ( 1.7) -( 1. 10)  that,  when  the  law  for  the 
steady-state  temperature  distribution  along  the  length  is  given,  the 
stability  of  the  flow  is  determined  by  three  dimensionless  criteria: 

1)7  —  the  specific  heat  ratio;  2)  Mo  —  the  inlet  Mach  number  of  the 
flow;  3)  p.  —  the  ratio  ( doubled  for  variant  2)  of  the  instantaneous 
heat  transfer  over  the  entire  length  of  the  channel  for  constant.  Qo 
(equal  to  the  heat  transfer  rate  in  the  inlet  cross  section)  to  the 
initial  heat  content  of  the  instantaneous  rate  of  flow  of  the  gas. 

The  parameter  p.  may  be  expressed  in  terms  of  the  ratio  of  the  initial 
and  final  temperatures: 

(3.») 

We  note  that  p.  characterizes  both  the  steady-state  flow  parameters 
and  the  heat  supply  rate  in  the  presence  of  oscillations  for  the  adopted 
quasi- steady  dependence  of  Q  on  p.  The  velocity  u0  enters  into  both 
p.  and  Mo.  Obviously,  when  uo  -*  0(p  °°  ,  Mo  -»  0)  the  steady-state 

regime  Is  not  possible. 

The  machine  computation  was  carried  out  for  several  values  of 
7  (l-£  y  4.  1.4)  ,  Ho  (0.01^  Mo^  O.l),  and  for  values  of  p.  from  zero 
to  70  and  above.  The  fundamental  characteristic  values  of  the  problem 
are  infinitely  numerous,  the  calculation  was  performed  for  the  first 
three  (k  =  1,  2,  and  3  in  (3.2))  . 

The  computations  showed  than  when  p.  is  increased  the  characteristic 
values  move  in  a  certain  way  in  the  xy  plane,  approaching  the  y-axis 
and  tending  toward  some  limiting  value  as  p.  -)  00 . 

The  least  damping  during  the  increase  in  p.  is  observed  in  the 
fundamental  note  of  the  oscillations  (k  =  l)  . 

For  variant  1  ( variable  temperature  gradient)  stability  boundaries 
in  the  p.7  plane  for  several  values  of  Mo  have  been  presented  in  Fig.  1. 
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Pig..  2.  Fig.  3. 


(Values  of  7  close' to  unity  are  possible  for  flows  with  dissociation 
or  ionization) .  The  frequency  of  the  oscillations  0  for  all  Mo  and  7 
on  the  stability  boundary  are  described  with  high  accuracy  by  Eq. 
(3-2)  when  VyCa 


A  curve  of  Cx\VlCa 


is  replaced  by  the  average  sound  velocity  C^: 


0 

V  5iv4 

is  given  in  Pig.  .2. 


(3-5) 
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Fig.  4. 


The  possibility  of  introducing  the  average  sound  velocity  indicates 
a  coincidence  between  the  velocity  of  the  propagation  of  the  oscilla¬ 
tions  and  the  velocity  of  sound  for  the  given  temperature  distribution 
The  velocity  of  the  propagation  may  be  different  when  other  distribu¬ 
tion  laws  apply  [3]. 

As  follows  from  Pig.  1,  the  flow  is  stabilized  during  the  increase 
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in  M0  and  7  and  the  decrease  in  p.  Presented  in  Fig.  3  are  the 
averaged  absolute  values  of  the  characteristic  functions  which  lie 
close  to  the  distribution  of  the  amplitudes  of  the  neutral  oscillations 
along  the  length  of  the  channel  {when  Mo  =  0.1,  7  =  1.2,  and  6(0)  =  l)  . 
In  contrast  to  isoentropic  oscillations,  for  which  tj  =  76,  the  ampli¬ 
tudes  of  the  density  (and  velocity)  oscillations  greatly  exceed  the 
amplitudes  of  the  pressure  oscillations.  Presented  in  Fig.  4  are  the 
calculated  real  and  imaginary  parts  of  the  density  oscillations 
6  -  +  162.  Present  in  the  curve  for  dj  are  the  oscillations  which 

are  superimposed  over  the  fundamental  oscillation.  These  oscilla¬ 
tions  have  a  shorter  wavelength  (by  about  Mo  times)  than  that  of  the 
fundamental  oscillation  and  correspond  to  isobaric  entropy  disturbances 
(these  oscillations  are  not  noticeable  in  the  curves  of  and  ^g)  . 

For  variant  2  ( constant  temperature  gradient  along  the  length) 
no  unstable  oscillations  were  detected  even  for  the  larger  values  of 
|i  (up  to  200)  .  This  attests  to  the  strong  influence  on  flow  stability 
of  the  steady-state  heat  supply  law.  The  flow  Is  also  stable  for 
q  =  0,  when  small  disturbances  do  not  alter  the  heat  supply  rate. 

Received  January  15,  19^2. 
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NON- STEADY  FLUTTER  OF  PLATES  AND  SLANTED  SHELLS 


IN  GAS  SHELLS 

V.  V.  Bolotin 


The  problem  of  the  unsteady  flutter  of  plates  -v 
and  slanted  shells  in  supersonic  gas  streams  is  , 
examined.  The  examination  Is  based  on  equations  \ 
which  were  obtained  in  previous  papers  [ 1  and  2] 
and  later  applied  to  a  number  of  steady  problems 
regarding  oscillations  and  thermal  buckling  [3-7]. 

The  present  article  deals  with  finding  non-steady 
solutions  for  the  simultaneous  variation  of  stream 
velocity,  ambient  medium  density,  plate  or  shell 
temperature,  etc.  For  this  an  approximate  method 
for  Investigation  of  non-steady  processes  In  linear  . 
systems  with  many  degrees  of  freedom  [8],  hereafter  \ 
referred  to  as  the  "method  of  single-frequency 
oscillations"  is  employed.  An  example  Is  given  of 
a  determination  of  the  amplitudes  of  non-steady 
flutter  for  a  flat  panel. 


1.  Initial  Equations.  By  the  method  described  in  previous 
papers  [i,  2,  and  6]  the  equations  for  the  oscillations  of  plates  and 
shells  lying  in  a  temperature  field  and  flowed  past  by  a  supersonic 
gas  stream  may  be  reduced  to  the  following  system  of  equations: 

ft  * 

ti  +  Sj£j  ~r  W;!{j  4-  0  2  ai*  £*  +  H  2  bjklk  =  (  1  •  l) 

*«=l 

(felt  (it  *  *  ■  »  CtiiCjt  (it  •  •  •  >  («t  gt  (/  “  1*2. .  .  .  .  il) 
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Here  £ 1  are  generalized  coordinates  ( the  dimensionless  coeffi- 
cients  in  the  expansion  of  the  additional  normal  deflection  of  the 
plate  or  slanted  shell  into  series  with  respect  to  the  n  first  modes 
of  the  natural  oscillations  in  a  vacuum;  g^  and  a>j  are  the  dimension¬ 
less  damping  characteristics  and  the  dimensionless  frequencies  of  the 
respective  oscillation  modes;  a^  and  are  certain  coefficients 
depending  on  the  nature  of  the  temperature  field  and  the  properties 
of  the  flow;  <f>j;  are  certain  linear  functions,  which  are  polynomials 
of  the  generalized  coordinates  and  generalized  velocities.  The  dots 
in  Eqs.  (l.l)  signify  differentiation  with  respect  to  dimensionless 
time  _t.  The  characteristic  panel  temperature  and  unperturbed  stream 
velocity  are  described  by  the  parameters  9  and  p  respectively.  We 
shall  choose  these  parameters  in  the  following  way  [J>  and  6]: 


0  =  -^  (1  +'W  (£)  . 


P  = 


(1-2) 


where  T  is  the  characteristic  temperature;  a  Is  the  thermal  expansion 
coefficient  of  the  panel  material;  v  is  the  Poisson  ratio;  E  is  the 
elasticity  modulus;  ji  is  the  characteristic  dimension  of  the  middle 
surface;  Ji  is  the  thickness  of  the  panel;  w  is  the  polytrope  index  of 
the  flow;  p.  is  the  pressure  of  the  unperturbed  flow;  and  M  is  the 
Mach  number  for  the  unperturbed  flow  (by  unperturbed  flow  we  mean  flow 
in  the  absence  of  flutter  or  thermal  buckling  of  the  paneling)  .  It 
is  possible  to  show  that  the  parameters  9  and  p  quite  completely 
describe  the  conditions  of  the  behavior  of  parceling  in  a  variable 
temperature,  velocity,  and  pressure  field.  Tt&e  values  of  9  and  p  at 
any  moment  of  time  yield  a  certain  image  point  on  the  9  and  p  plane. 

The  linear  system  corresponding  to  (l.l) 

n  % 

w j 4-  £jt j  +  ®j*£>  +  0  2  aj*t*  + I* 2  =  0  —  . B)  (l,j) 
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has,  apart  from  the  trivial  solution  £  1  *  »  .  ..^£n  =  0,  a 

solution  of  the  type 

£'”E'e“  (1.4) 

where  and  a  are  certain  constants.  Substituting  (1.4)  into  Eq. 
(1.3)  we  obtain  the  characteristic' equation 

|  (a*  +  g,s  +  u/)  6  jj,  +  0  ajk  +  |  =  0  (1.5) 

(6^  is  the  Kronecker  symbol).  The  trivial  solution  is  stable  as 
long  as  all  Re  a  <  Oj  it  becomes  unstable  if  Just  one  of  the  roots 
of  Eq.  (1.5)  passes  over  into  the  right  half-plane. 


Pig.  1. 


Hereafter  it  will  be  necessary  to  distinguish  between  two  types 
of  transitions  to  instability.  First,  the  transition  may  take  place 
through  a  point  on  the  imaginary  axis  not  coincident  with  the  origin 
(Fig.  l).  In  this  case  Eq.  (1.5)  has  on  the  stability  boundary  two 
purely  imaginary  roots  a  =  +  ico#,  where  oj*  is  the  flutter  frequency. 
On  the  p.0  plane  the  set  of  points  corresponding  to  this  type  of 
instability  form  a  line  f(li*,9*)=  0,  which  we  shall  call  the  line 
of  dynamic  instability  or  the  flutter  line.  Second,  a  transition 
through  the  origin  is  possible  (Fig,  2).  In  this  case  Eq.  (1.5)  has 
a  root  a*0,  and  the  stability  loss  is  of  a  non-oscillatory  nature. 


The  line  s(p#,  0.J  *  0  on  which  one  of  the  roots  of  Eq.  (1.5)  has 
the  root  0  »  0,  while  the  remaining  roots  lie  in  the  left  half-plane 
will  be  called  the  line  of  static  instability. 


The  arrangement  of  these  lines  for  the  case  of  a  flat  rectangu¬ 
lar  plante  [3]  is  shown  in  Pig.  3. 

When  the  transition  takes  place  through  the  line  AF,  flutter  en¬ 
sues,  and  when  it  takes  place  through  the  line  DF,  a  thermal  buckling 
of  the  plate  occurs.  Corresponding  to  a  variation  in  flight  condi¬ 
tions  are  fairly  complex  trajectories,  which  the  image  point  describes 
on  the  |X0  plane. 

2.  Use  of  the  Method  of  Single-Frequency  Solutions  for  Analysis 
of  the  Behavior  of  a  System  Close  to  the  Flutter  Line.  Let  us  assume 
that  the  image  point  intercepts  the  flutter  line  f(n*,0.J  =  0  at  some 
point  N  (Fig.  3)  •  Let  be  the  frequency  of  the  oscillations  at  the 
flutter  boundary.  If  the  boundary  is  crossed  slowly  enough  and  if 
the  image  point  does  not  go  too  far  from  the  boundary,  then  it  is  to 
be  expected  that  the  oscillations  of  the  system  will  be  close  to 
periodic  oscillations  with  frequency  o>#. 

We  shall  use  this  idea  as  the  basis  for  an  approximate  solution 
of  system  (l.l)  for  variable  p.  and  0. 

After  writing  system  (l.l)  in  the  form 
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(2.1) 


ii  4  Sjlj  +  «i*£j  +  6  (t)  2  ajkU  4-  M(T)  2  *;*£*  = 

*=•  *-i 

=  ®j  (Cl.  £*••••  £*.  Cl.  Cl. - in,  P  (t),  0  (t)J  (/  =  1,2 . n) 


we  shall  assume  that  the  parameters  0  and  p  are  functions  of  the 
"slow"  time  t '»  et  ( e  «  l)  .  Furthermore,  we  shall  assume  that  the 
"mistunlngs"  p-p*  and  0-0*  are  sufficiently  small.  For  the  time  being 
we  shall  make  no  assumption  concerning  the  smallness  of  the  damping 
characteristics  g^j  as  we  shall  see  below,  this  is  because  of  ce_  tain 
peculiar  features  of  the  problem  under  consideration.  As  before,  the 
dots  in  Eqs.  (2.l)  signify  differentiation  with  respect  to  "fast" 
time  _t. 

Let  us  transform  Eq.  (2.l)  in  the  following  way: 

*  * 

t>  -!-  Pit)  4-  “ftj  -T-  0.  2  +  H.  2  fy*C*  «■ 

-  •  *-i 

=  eVyl?,,  £s . £„  Cj,  . . j„,  n(r),  9(t)J  {,  =  1^ . B) 

eM^  =  (Vj  !£„  . U  £s,  C, . U  |i  (t),  0  (t)J  - 


(2.2) 


where 


-  (8  -  0.)  2  «j&  -  (I*  -  P.)  2 

*=i  *-i 


(2.3) 


We  introduce  here  the  small  parameter  e  in  accordance  with  the 
above  made  assumptions.  It  may  be  assumed  that  the  point  N  "wanders" 
along  the  flutter  line:  0*  =  0*(t),  p*  =  p*  (t).  For  example,  it 
Is  advantageous  to  choose  It  so  that  the  overall  mistuning  [(p  —  p*)  2+ 

+(0  —  0*)  2]  1/^z  be  minimized  (Fig.  4),  etc.  But  here  for  simplicity 
we  shall  assume  that  p*  =  const  and  0*  *  const. 
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Let  us  consider  the  generating  system  corresponding  to  (2.2) 


tj  -+-  gj  t)  -t-  ®i*  t j  -r  0.  V  ajgt  +  j».  V  />,»£*  =  0  </  -=  1,2 . ») 

*  - 1  *1 

We  note  the  following  characteristics  of  system  (2.4): 

1)  system  (2.4)  admits  of  two  linearly  independent  periodic 
solutions  with  frequency  co*.  In  fact,  Eq.  (1.5)  has  two  purely  . 
imaginary  roots  a  =  +  to#  on  the  flutter  boundary.  Consequently,  it 
is  possible  to  construct  a  solution  to  system  (2.4)  as  a  function  of 
two  real  constants. 

2)  The  only  steady-state  solution  to  system  (2.4)  is  the  solu¬ 
tion  c  1  -  Ca  •  ••  •  *  Cn  -  0. 

3)  The  frequency  id*  is  not  multiple,  i.e.,  at  the  point  N 
(Fig.  l)  only  one  of  the  characteristic  indices  emerges  onto  the 
right  half-plane.  We  note  that  if  the  system 

\  n  n 

T  t  8.  2!  +  F.  i 

*-l  »- 1 

were  taken  as  the  generating  system  (which  is  natural  in  view  of  the 
relative  slightness  of  the  damping)  ,  this  condition  would  not  be 
fulfilled.  This  peculiarity,  which  is  illustrated  in  Fig.  5,  wad 
first  noted  in  one  of  the  author's  previous  papers  [2].  The  condition 
of  the  absence  of  multiple  frequencies  is  a  portion  of  the  stronger 
condition  of  the  absence  of  "inherent  resonances".  But  if  the  degree 
of  instability  is  equal  to  two;  i.e.,  if  all  the  characteristic  indices 
except  two  lie  in  the  left  half-plane,  the  condition  that  has  been 
formulated  is  sufficient. 

It  is  easy  to  see  that  the  theory  of  "single-frequency 
oscillations"  in  systems  with  many  degrees  of  freedom  developed  by 
N.  M.  Krylov,  N.  N.  Bogolyubov,  and  Yu.  A.  Mitropol '.skiy  [8]  is  entirely 
applicable  to  system  (2.2).  The  generating  solutions  are  of  the  form 
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W"  = 


(2.5) 


(the  superscript  asterisk  denotes  conversion  to  the  complex  conjugate). 
The  complex  constants  are  found  as  a  non-trivial  solution  to  the 
algebraic  equations 

»  n 

<“i*  +  -  <■>.*)  li  +  0.  23  -f  K  2  =  0  (/  »  1 .2  ..  „)  (  2  • 6) 

We  shall  seek  a  solution  to  the  non-linear  system  (2.2)  close 
to  the  solution  to  (2.5)  in  the  form 

Cj  +  “li'e-**  +  eu/u  (a,x p)  +  e2Uj<*>  (a,  \p)  +  . . .  (2.7) 

Here  a  and  q>  =  f  —  oj^t  are  the  slowly  varying  amplitude  and 
phase  (real  functions  of  the  "slovi"  time  t)  ,  and  Uj^  Uj^  , . . .  are 
the  amplitudes  of  a.  The  functions  a(r)  and  f(i)  are  determined  from 
the  equations 

e-4j (a) +  r.M, («)  +  ... 

5T=  *.  +  eIi,  (•>)  +  t2Bz  (a)  +  . . .  (  2-8) 


i 

I  mt 

•e» 

_ _ i 

0 

n 

I  ^ 

, 

Fig. 
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It  is  necessary  to  set  up  an  algorithm  for  finding 

u^  ^  ,  Uj^  ^  ,  . . . ,  Ai,  A2, ...,  Bi,  Bfi,  . ..,  such  that  the  expressions 

(2.7)  satisfy  Eqs.  (2.2)  whenever  a  and  ^  satisfy  Eqs.  (2.8). 

Hereafter  we  shall  confine  ourselves  to  the  first  approximation 
The  equations  for  the  first  approximation  can 
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be  derived  from  (2.2)  by  a  somewhat  simpler,  though  more  formal, 
method  than  was  done  in  Bogolyubov's  paper  [8].  In  the  present 
article  we  shall  choose  this  simpler  method. 

3.  Derivation  of  the  Equations  for  Amplitude  and  Phase  in  Non- 
Steady  Flutter.  In  contrast  to  (2.7)  the  solution  to  (2.2)  will  be 
sought  in  the  form 

Ej  = 

where  £  and  f  are  real  functions,  and 

%  =  Ej  +  (a.  V)  +  («,*)  +  ••• 

%  =  £;*  +  *V’  (a.  ^)+  (a,  4)  +  . . . 

If  the  desired  solution  is  close  to  a  harmonic  with  frequency 

then  expressions  (3-2)  may  obviously  be  treated  as  functions  of 

"slow"  time  t.  Taking  this  into  account  we  may  write  down  the  following 

approximate  relationships: 

Ei  =  la  +  ‘‘a  (<■>.  +  4)1  +  |a  —  ia  («.  +  q>)]  +  . . . 

Ej  =  [2i'a <d.  —  a  (a>.5  +  2o).  <p)J  +  (3.3) 

•f  (—  2ia  o».  —  a  (u>.*  +  2o>.  «p)]  rjja  e~**  -f-  . . . 

The  dots  signify  terms  containing  c2,  e3,  etc.  Expressions 
(3.3)  are  substituted  into  the  left  side  of  Eqs.  (2.2).  The 
expressions 

c,  =  +  ag;  <r^  +  ...  &  =  iam^*  —  law.gj*  tr»  +  ...  (5.4) 

are  substituted  into  the  right  side,  since  it  contains  a  factor  e. 
Expanding  the  right  side  into  a  Fourier  series 

V}  =  Y/'V*  +  VjM  er**+  ...  . 

are  equating  coefficients  exp  If,  we  obtain  the  system  of  equations 

(<Dj‘  -  <0  arjii  +  |«».  emu  +  0  a  +M.a  J  =  (5.5) 

■  ■  ■  -  .  *”* 

=  —  2i0).  ali  +  2w.«f«5;  —  g^  —  +  tp^ 


(5.1) 

(3.2) 
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Equating  coefficients  of  exp(— 1 f)  obviously  leads  to  an  equivalent 
set  of  equations.  System  (3-5)  may  be  formally  treated  as  a  system 
of  linear  algebraic  equations  relative  to  tj^.  By  virtue  of  the 
choice  of  the  generating  solution 

A  =  I  (O)/  —  0).*  +  igpj  +  O.flj*  +  n.6*  I  =  0  (3.6) 


In  order  for  system  (3**0  to  admit  of  a  solution  it  is  necessary 
and  sufficient  that 

n 

2  (—  2i<0.  ali  -f  Soj.qwg, —  +  'Ff>)Xj*=0  (3.7) 


Here  Xj  is  a  solution  to  a  system  conjugate  with  respect  to 
system  (2.6).  Let  us  denote 

s  vx  s  i&x; 

^ - =  G+iH,  t± - =  g 

S  2  Mi 

>=i  j-i 

where  G  and  H  are  real  functions  of  a.  Condition  (3.7)  assumes  the 
form 

—  (2i®.  +  g)  o  +[(2®.  -  ig)ay  =  G  +  iH  (3-8) 


By  separating  the  real  and  imaginary  parts,  we  obtain  from  here 


2ci),//  (a)  +  gG  (a) 
4<1>.*  -j-  g*  * 


= 


2<a,G  (a)  —  gH  (a) 

a  (4u.*  +  *») 


(3.9) 


For  small  £  we  arrive  at  the  ordinary  equations  of  the  Krylov- 
Bogolyubov  method 


9  =  «.  + 


2aom 


(3.10) 


We  obtain  the  steady  solutions  by  setting  a  *  0,  9  *  const.  The 
equations  for  this  case  obviously  coincide  with  the  equations  obtained 
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by  the  small  parameter  method  [2]. 

Instead  of  using  the  conjugate  system  concept,  it  is  possible 
to  proceed  in  the  following  manner.  Let  us  consider  the  determinant 
of  (3.6)  denoted  by  A.  The  generating  system  admits  of  a  unique 
periodic  solution  with  a  simple  exponent  equal  to  ia>#;  therefore 
the  rank  of  the  determinant  is  obviously  n-1.  In  order  for  system. 
(3-4)  to  be  solvable  with  respect  to  aq^,  ^  necessary  and  suffi¬ 
cient  that  the  determinant  Aq,  which  is  obtained  from  A  by  replacing 
one  of  the  columns  by  the  column  of  the  right  side  of 

Vj  =  —  (2uo.  a  —  2u.(f  a  +  gp  +  igjipa)  ls  +  ¥/*?  (  3 . 1 1) 

be  equal  to  0.  After  separating  the  real  and  imaginary  parts  we 
obtain  the  equations 

Be  A0  =  A,  (o',  9,  a,  9,  n,  0)  =  0  (  5 . 12) 

Im  A0  =  A.  (a,  9,  a,  9,  ji,  6)  =  0 

equivalent  to  Eqs.  ( 3 - 9) - 


4.  Example.  Flat  Unheated  Panel.  Let  us  assume  that  a  flat 
rectangular  plate  with  sides  a  and  b  and  constant  thickness  h 
supported  along  the  perimeter  is  flowed  past  by  a  supersonic  flow 
in  the  direction  of  one  of  the  sides  ( Fig.  6)  .  We  shall  proceed  from 
the  nonlinear  Karman  equations,  neglecting,  as  is  usually  done,  the 
tangential  components  of  the  inertial  forces  and  considering  the 
perimeter  of  the  panel  indc._  ormable.  The  initial  stresses  in  the 
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middle  surface  will  be  assumed  equal  to  zero.  The  excess  aerodynamic 
pressure  will  be  determined  from  the  linearized  quasi- steady  formula 
of  '’piston”  theory.  The  deflection  w(x,  y,  t)  will  be  approximated 
with  the  aid  of  two  modes  of  oscillations 

[.jut  _  .  iT*!  .  xy 

£1  <0  7-  +  £*  (0  »n  —  J  *  *»“  -jf- 

Under  the  assumptions  that  have  been  made,  system  (l.l)  has  the 

form 

Ci  +  +  “i*»i  —  *3*  (£1.  £1).  f+  f»£i  +  “t*£i  +  "J"  |*£i  =  (£1.  £>)  (4.l) 


Here 


-[ 


1  TO*  — 


“i«=  1  +  »*,  wi“4  +  m'  (mms  — J 
®i  =  ~  £1  (<  1  i£i*  +  <‘h£i5).  =  —  £1  (tn£i*  +  CbJj’) 

2  (1  -f  2vm*  J-  m4 


-V* 


fn  =  (1  +  m*)  + 


lc-  T..  .  .  ,  2(16  +  8vm«  +  m41  _ 

Jo,  Cn~pC-f'n<-i - j  _vi - J  .S 

’  > 


81m4  m4 

(TV  4m-)«  +  (9  +  4to»)»  + 


2(4  +  5vm»+m*)l 
l-v» 


S  =  T(>  —  v*) 

For  the  dimensionless  time  t  we  shall  take  the  parameter 


(4.2) 

(4.3) 

(4.4) 


n*  [  D  \  » 
=  «*  (p»a) 


Here  ti  is  natural  time;  D  is  the  cylindrical  rigidity;  and 
Po,  the  density  of  the  panel  material.  In  subsequent  calculations 
we  shall  assume  gx  =  ge  =  g. 

It  is  not  difficult  to  deduce  the  following  formulas  for  the 
critical  parameter  p.*  and  the  frequency  of  the  oscillations  at  the 
boundary  of  the  flutter  region  ®#: 


j*.  =  l(®i*  —  <oi*J'  4-  2 g*  («,*  +  «*») jv',  ».  =  )  ( ^  •  5) 

and  also  expressions  for  the  constants  £1  and  £2.  Since  these  constants 
are  related  by  the  formula 

y  J*.  5i  +  {«•'—«.*+  igu, )  «=  0 
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we  may  take  for  them  the  expressions 

{.  =  4-*  l»  =  —  y  «-,T 


Here  7  is  the  phase  shift  between  £1  and  — |a 

ttnT (4.6) 

For  such  a  choice  of  the  parameters  and  £z  an  amplitude  a  =  1 
corresponds  to  a  flexure  along  one  halfwave  in  a  direction  parallel 
to  the  flow  with  a  maximum  deflection  equal  to  the  thickness  of  the 
plate  h.  The  amplitude  of  the  other  component  in  this  case  will  also 
be  equal  to  h  (Fig.  7)  . 


Fig.  7. 


Fig.  8. 


Proceeding  to  construct  Eq.  (3.9) ,  we  shall  first  write  out 
the  determinant  in  (3-6) .  In  the  notations  of  (4.6)  it  has  the  form 


(4.7) 


Consequently,  the  condition  Aq  *  0  obtained  from  (4.7)  by  re¬ 
placing  one  of  the  columns  by  a  column  made  up  of  the  functions  (3-10) 


will  be 


i>iT  +  v, = 0 


(-4.8) 
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Let  us  calculate  the  functions  Vx  and  V2.  By  noting  that,  when 
Eqs.  (3-4)  and  (3-5)  are  taken  Into  account, 

:.*=  4  •»(«<*+ 

W  =  -  T  “*  I'1  <*~T)  +  +  •  •  • 

=  -  T  «*  +  *<Y)  *'*  +  (2*W  +  e-*T) e-*|  + . . . 

=  T •'  1(2  +  * - *Y)  ***  +  (2  +  e**T)  «*]  +  ... 

(dots  denote  the  terms  containing  the  harmonics),  we  find  from  ( 3 . 1 l) 
and  (4.2) 

2  (u  —  u. )  a  ,v  3a*  I  2  +  e— **T  V 

2V,  =  —  (2jci),o  —  2u, tfa  +  ga  +  ig( fa)  — - 3 -  e  T  —  -y'  (eu  + jj - r„ J 


21’,  =  (2ib>.a  —  2ai.(ja  -j-  go  +  '*9°) «  ,Y— 


-iT  2  (n  —  H.)  a 


3a*  1 2c~ <T  +  c*1 
‘  4  \  3 


ea  +c»c- 


Substitution  into  Eq.  (4.8)  yields 

-(2i*».  +  *)«-$- <2».—i s’n  y  -  [y  (,i  -ji.)a  - 
— 4*  o’  (<•«  -f  fa  —  Cu  —  Cu)  cos  T  +  -|-  o’l  (r„  -f  e»  +  — y*”)  sinr] 

Combining  this  equation  with  (3-8),  we  find  that 

C<“>  =  T“*  (cu  +  c0  +  ^5i±fE) 

II  (a)  =  _  1  T-i /  3  l  (4.9) 

2  sin  T  L  3  H.)  a  —  J  o’  (rn  +  c»  —  e,,  —  e„)  cos  jl 

To  obtain  the  steady-state  amplitude  a  =  ao  =  const  we  have  the 


equation  2co  H(ao)  +  gG(a0)  =  0.  Hence  we  arrive  at  the  formula 


*  I*- M% 

-“tv ~ur) 


(4.10) 


where  c0  =  c2X  +  c22  —  Cn  —  cx2,  and  the  coefficient  k  is  defined 


by  the  formula 


,  _  (j  ,  f  tarf  fit  +  cn  +  (fii  +  f»)  / 3  ' 
V 1  "**  2a>.  fn  +  fn  —  cu  —  e„  , 


-  (fii  +  f»)  /  3  \ 

—  CU  —  Cn  /  ‘ 


(4.11) 


As  is  apparent  from  (4.6) ,  tan  7  >  0,  tan  7  being  about  equal 
to  £.  If 


f'^ui'-tti* 


(4.12) 
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we  may  assume  that  k  *  1,  and  Eq.  (4.10)  converts  to  the  well  known 
formula  encountered  In  the  literature  [1,  2,  6J. 


Pig.  9.  Pig.  10. 


Assume  condition  (4.12)  is  fulfilled.  The  fii*st  equation  in 
(3-9)  takes  the  form 

-TsSirr^-*)  (4.13) 

Shown  in  Fig.  7  are  the  three  ideal  programs  for  the  variation 
of  the  parameter  p.  in  time.  The  results  of  numerical  integration 
of  Eq.  (4.13)  with  the  initial  condition  a( t*)  *  0.01  are  given  in 
Fig.  8-10  ( t*  is  the  moment  of  time  corresponding  to  the  first  inter 
section  of  the  flutter  line) .  It  is  assumed  that  the  plate  is  a 
square  (m  =  l)  ,  while  the  Poisson  constant  v  =  0.3.  Under  these 
assumptions 

***•  ~ 15,75'  “.  ~  3.808,  c,  =  35,438. 

The  curves  of  a  =  ao( t) ,  which  correspond  to  the  instantaneous 
establishment  of  steady  solution  are  plotted  by  the  broken  lines. 
Calculations  carried  out  for  the  cases  sin  7  *  0.01,  0.05,  and  0.10  . 
demonstrate  the  strong  effect  of  damping  on  the  amplitude  of  the 
non-steady  flutter. 
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